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Abstract 

In this paper we construct a covariant light-front model of heavy mesons 
within the framework of heavy quark effective theory (HQET). The covariant 
model consists of the light-front heavy meson bound states constructed in 
the heavy quark limit with heavy quark symmetry and heavy quark effective 
theory as its basis. Within this model, the Isgur-Wise function and decay 
constants in the infinite quark mass limit can be evaluated in a very simple 
and the most general way. The results are ensured to be consistent with heavy 
quark symmetry. From the heavy-quark-limit bound states, we can further 
develop a systematic approach to calculate 1 / tuq corrections from the 1 / mq 
expansion of QCD. This covariant model can serve as a quasi-first-principles 
description of heavy meson dynamics, namely, a phenomenological covariant 
bound state in the heavy quark limit which is consistent with heavy quark 
symmetry, plus a reliable first-principles computation of the 1 /tbq corrections 
in HQET in terms of the 1/mn expansion of the fundamental QCD theory. 



PACS numbers: 12.39.Hg, 12.39.Ki, 12.60.Rc, 14.40.-n 
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I. INTRODUCTION 



In the past decade, the most significant progress made in the QCD description of hadronic 
physics is perhaps in the avenue of heavy quark dynamics. The analysis of heavy hadron 
structures has been tremendously simplified by heavy quark symmetry (HQS) proposed by 
Isgur and Wise IJ and the heavy quark effective theory (HQET) developed from QCD in 
terms of I/tuq expansion 0. HQS predicts that all the heavy to heavy mesonic decay 
form factors in the infinite heavy quark mass limit are reduced to a single universal Isgur- 
Wise function, while HQET provides a systematic framework for studying of the symmetry- 
breaking I/frig corrections (for a review, see ||). Applying HQS and HQET to heavy 
baryons leads to similar simplifications to heavy baryon structures |3J]. Also, combining 
with chiral perturbation theory, one can construct a chiral Lagrangian for the low-energy 
interaction of heavy hadrons with Goldstone bosons |J. Moreover, in terms of heavy quark 
expansion, HQET offers a new framework to systematically study the inclusive decays of 
heavy mesons ||. 

However, the general properties of heavy hadrons, namely their decay constants, tran- 
sition form factors and structure functions etc., are still uncalculable within QCD, even in 
the infinite quark mass limit with the utilization of HQS and HQET. Indeed, HQS itself is 
unable to determine the explicit form of the Isgur- Wise function and it does not and can- 
not provide a direct description of nonperturbative QCD dynamics to other heavy hadron 
properties. The lattice QCD simulations J7J permit a nonperturbative approach to solve 
the nonperturbative QCD problem, but so far a direct calculation with the 6-quark is still 
not possible due to the difficulty of placing heavy particles on the lattice. Very recently an 
alternative attempt in this regard has been carried out from light-front QCD || but further 
investigations are needed to obtain reliable results. Hence, although HQS and HQET have 
simplified very much the heavy quark dynamics, a complete first-principles QCD description 
of heavy hadrons is still lacking due to the unknown nonperturbative QCD dynamics. 

In most studies, the heavy hadron decay constants and transition form factors are usually 
evaluated using various phenomenological models, such as the constituent quark model 0, 
the MIT bag model [ IU| , QCD sum rule [IT], and the light-front quark model Ji"2|-|2"2]] . The 



nonrelativistic constituent quark model and the MIT bag model have been widely used in the 
phenomenological discussion of hadronic structures. However, the application of these two 
models is trustworthy only for the processes with small momentum transfer. The problem 
may be overcome by the light-front quark model (LFQM) [|12|] which has been considered as 
one of the best effective relativistic quark models in the description of the exclusive heavy 
hadron decays [|13H 22|1. Its simple expression, relativistic structure and predictive power 
have made wide applications of the LFQM in exploring and predicting the intrinsic heavy 
hadron dynamics. However, theoretical justification of such a model within the framework 
of QCD has not been addressed in the literature. 

In fact, the phenomenological LFQM is also largely limited in its applications, due mainly 
to the quark model assumptions and its non-manifestation of covariance as a relativistic 
description. For example, it has been emphasized in the literature that various hadronic 
form factors calculated in the LFQM should be extracted only from the plus component 
of the corresponding currents because of the dynamical dependence of other components. 
Also, the LFQM calculations for decay processes are restricted only with the zero momentum 
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transfer. As a result, one often picks up a specific Lorentz frame to perform such light-front 
calculations. These computations can lead to the correct results only if the covariance is 
precisely preserved within the LFQM formulation. We have recently attempted to calculate 
the decay form factors over the whole range of momentum transfer, and realized that the so- 
called Z-diagram (a higher Fock state contribution beyond the LFQM) must be incorporated 



in order to maintain the covariance ||20|| . However, with such an extension, the simplicity 
of the LFQM is lost. Without including the Z-diagram contribution, the heavy hadronic 
form factors calculated directly in the LFQM over the whole range of momentum transfer 
will explicitly break relativistic covariance JH],^. Thus, the current approach to extract 
hadronic form factors in the LFQM formulation by choosing a specific Lorentz frame and 
then calculating a particular component (the plus component) of the associated current 
matrix element may not be unique and may cause some inconsistent and even misleading 
results. 

To resolve the above-mentioned problems in the LFQM and to improve the current 
understanding of the QCD analysis of heavy hadrons, we construct in this paper a covariant 
light-front model. This model consists of a heavy meson bound state in the heavy quark 
limit, which is phenomenological in nature at the present time but is fully consistent with 
HQS, plus a reliable approach from this bound state to systematically calculate the l/m,Q 
corrections within HQET in terms of the 1/rriQ expansion of the fundamental QCD theory. 
We expect that this covariant model fully based on HQS and HQET can serve as a partially 
phenomenological (the bound states in the HQS limit) and partially fundamental (the I/tuq 
corrections in QCD) description of heavy mesons. 

The paper is organized as follows: We will provide in Sec. II a general construction of 
light-front bound states within HQET. To have a consistent light-front description of heavy 
hadrons, we begin with an explicit boost covariant expression of meson light-front bound 
states which can be further expressed in a covariant form in the heavy quark limit. We 
show from the light-front bound state equation that such bound states obey HQS. With 
the covariant bound states in the heavy quark limit, we can further develop a systematical 
approach to calculate the 1 / vtlq corrections to the heavy meson energy splitting, decay con- 
stants and transition form factors within HQET. This will provide a reliable first-principles 
QCD analysis of the l/m^ contributions to heavy hadron dynamics. 

For a consistent check, we evaluate in the Sec. Ill the Isgur-Wise function £(v ■ v') in a 
complete covariant way so that no specific Lorentz frame has to be pre-fixed, and it is not 
necessary to only calculate the plus component of the current in matrix elements. We can 
straightforwardly show that in the covariant model the Isgur-Wise function is renormalized 
to unity at the zero-recoil point without using the detailed form of the wavefunction, as 
predicted by HQS. The general result of the Isgur-Wise function is found as a universal 
function of heavy meson decays. We also evaluate the heavy hadron decay constants Fp 
and Fy in the heavy quark limit in a covariant way. Again, a consistent result with HQS is 
obtained. 

In Sec. IV we further examine the covariant requirement for the light-front wave functions 
and calculate the Isgur-Wise function and decay constants. In Sec. V we evaluate for the 
first time the heavy mesons mass splitting in the covariant model. 

We can see that the covariant light-front model has overcome the theoretical difficulties 
encountered in the conventional light-front formulation for heavy hadrons, removed the 
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ambiguities in previous calculations, largely simplified the procedure of the non-covariant 
light-front formulation currently used in the literature, and further provided a first-principles 
QCD analysis of the I/tuq corrections within HQET. 



II. A CO VARIANT LIGHT-FRONT MODEL IN HQET 

The general form of the phenomenological light-front hadronic bound states that are 
currently employed in the literature has a similar structure as those in the constituent quark 
model and preserves a fully relativistic description. Hence these states are often called the 
relativistic quark model or the light-front quark model. However, almost all the previous 
investigations have not carefully treated the Lorentz structure of these bound states and 
paid enough attention to the consistency with heavy quark symmetry and heavy quark 
effective theory. As it is well known, a great advantage of the light-front bound states, 
which offer a fully relativistic description, is the simply manifested boost transformation in 
these bound states. To correctly describe the heavy quark dynamics, one should also ensure 
that the light-front description can reproduce all the results associated with HQS. Also, to 
accurately extract the heavy hadron form factors, the covariant problem of the light-front 
description of heavy hadron bound states must be considered. Furthermore, incorporating 
HQET with the phenomenological bound states in the heavy quark limit to analyze the 
1/rriQ correlations is strongly desired. These are the main issues we shall address in the 
development of our covariant model. 



A. Boost invariant form 

To make the boost symmetry manifest, we begin with the mesonic bound states that are 
written in a form of exhibiting explicitly the boost covariance, 



|P + , P ± , S, S z ) = J [d 3 Pl ][d 3 p 2 ]2(2n) 3 P + 5 3 (P - Pl - p 2 ) 



x 



^ SSz {x,K ± ,X 1 ,X 2 )\p 1 ,X 1 ;p 2 ,X 2 ) , (1) 

Ai,A2 

where we have used a boost invariant measure and the corresponding boost invariant light- 
front 5-f unction: 

[d 3 p] = > 2 (2^) V^ 3 (P - p') = 2(2tt) V<KP + - p' + )S 2 (p ± - p' x ). (2) 

In Eq. (|l|), pi,p 2 and Ai,A 2 are the momenta and helicities of the valence quark and an- 
tiquark, respectively, carried in mesons, and ^ SSz is the amplitude of the corresponding 
valence qq sector (and is simply called the light-front wave function) which depends only on 
the longitudinal momentum fraction x of the valence antiquark and its relative transverse 
momentum k± with respect to the center-of-mass frame of the hadron, 

X = , K L = p 2 ± - xP± . (3) 
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It is obvious that (x, k±) is boost invariant. In general, ty SSz can be expressed as 

tf 55 *(x, k ±} Ai, A 2 ) = FR ss *(x, K± , Ai, A 2 )$ s& (x, 4) , (4) 

with the flavor part of the wave function which is the same as that in the constituent quark 
model, and R and $ the spin and space parts, respectively, which depend on the quark and 
gluon dynamics. Note that R is also momentum dependent due to the relativistic feature 
of quarks. In phenomenological calculations, one usually ignores the dynamical dependence 
of the light-front spin so that R can be approximately expressed by the so-called Melosh 
matrix [23]. The valence quark Fock space \pi, Xi',P2, A 2 ) is given by 

|Pi,Ai;j?2,A 2 > =6 t (pi,Ai)d t (p 2 ,A a )iO), (5) 

with 

{b(p, A), b\p\ A')} = {dip, A), dV, A')} = 2(2vr) 3 p + 5 3 (p - p')8 xx , , (6) 

which is also boost invariant. 

In principle, the above light-front wave functions can be dynamically determined from 
the truncated light-front bound state equation J25| (which is similar to the Bethe-Salpeter 
equation in the usual equal-time frame): 

Ml-M^y^(x,K ± ,X l )= J ^^E^//(x,^,A j; x',4,A0^(x',^,A:), (7) 



where Mh is the hadron mass, Mq = ( P i + p 2 ) 2 = K \V^ 1 + K± ^ m " 2 is the so-called (boost) 
invariant mass obtained from the free energies of the constituents in hadrons |25[, and V e ff 
denotes an effective two-body interaction kernel. Note that by "truncating" the Fock space 
to only the valence quark states, the dominant contribution of higher Fock space to the 
bound states is described effectively by V e //- Such an effective two-body interaction on the 
light-front can be derived directly from QCD, for example, by the use of the light-front 
similarity renormalization group approach developed recently |26| . 



B. Heavy quark limit 

In this paper, we will focus on pseudoscalar and vector heavy mesons (corresponding 
to S — and 1, respectively) in the heavy quark limit. Let P = M^v in Eq. ([!]), where 

(v 2 = 1) is the velocity of the heavy meson and Mu its mass, and rescale the bound 
states (|l|) for heavy mesons by \P + , P±, S, S z ) = \/Mh\H(v, S, S z )). Then in the heavy 
quark limit, namely mq — > oo, Eq. ([I]) can be recast into 

\H(v,S,S z ))= J[d 3 k}[d 3 p q }2(2n) 3 v + 6 3 (Kv-k-p q ) ]T R SS *(X, K± , X Q , X q ) 

x<!> ss *(X,Kl)bl(k,X Q )dl(p q ,X q )\0), (8) 

where A = Mh — vtiq is the so-called residual center mass of heavy mesons, and = 
Pq — rnqv^ the residual momentum of the heavy quark with heavy quark mass rriQ. Here 
we have also defined the new light-front relative momentum: 
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< X = /v + < oo , — oo < k± — p q ± — Xv± < oo . (9) 

Note that we use the variable X (instead of the usual longitudinal momentum fraction 
x = p q /P + , where P^ = M^v^, or the residual longitudinal momentum fraction y = p+/K + , 
where = At> M ||) because it is the appropriate longitudinal variable appearing in the 
light-front heavy meson bound state equation in the infinite quark mass limit, see Eq. (|T8|) . 
The variable X was first introduced in [27] as X = Mhx, and it is also related to the residual 
momentum fraction y by X = Ay. The operator b\{k, Xq) creates a heavy quark with the 
residual momentum k^ and helicity Xq, 

{b v (k, X Q ), b\,{k', X' Q )} = 2(2n) 3 v + 5 vv ,5 3 (k - k')5 XQ y Q , (10) 

and [d 3 k] ee 

The spin part in the above bound state can be simply constructed in terms of the Lorentz 
structure of the pseudoscalcar and vector mesons: h v (x)i'j 5 q(x) and /i^(x)7 M g(x) whose 
momentum structures are R SSz proportional to u{v , Xo)irfv{p q , X q ) and —u(v, Xo)^v{p q , X q ), 
respectively, where the imaginary number i and the minus sign are introduced by a suitable 
choice of the phase factor to the pseudoscalar and vector meson bound states, u{v , Xq) and 
v(p q , X q ) are spinors for the heavy quark and light antiquark, 

J2u(v,X)u(v, X) = 1 + i , ^v{p q ,X)v{j) q ,X) = i> q -m q , (11) 

A A 

and e M is the polarization vector of the vector meson. Then the normalized R SSz is given by 

\ \\l — u(v, Xo)il b v(p a , X a ) for S=0, 
R SS *(X, k ± , X Q , X q ) = { 2 V 7 g+ 7 9 _ Q) K q q ' (12) 

l-2\/^+^ u(v,X Q )/v(p q ,X q ) for S=l. 

This result is indeed the heavy quark limit of the Melosh transformation given in Eq.(|4]). 

The space part $ 5Sz (X, k\) (called the light-front wave function) in Eq.(|8|) is the heavy 
quark limit of $ 5Sz (x, k±) in (^) rescaled by a factor \/Mh- The normalization of the heavy 
meson bound states in the heavy quark limit is then given by 

(H{v', S', S' z )\H(v, S, S z )) = 2(2nfv + 5 3 (Av' - Av)5 S s>8 Sz s>, (13) 
which leads to the following wave function normalization condition: 

dXd 2 n ± SSz 



2(2vr) 3 X 



\^(X,k 2 ± )\ 2 = 1. (14) 



C. Consistency with HQS 

In principle the heavy quark dynamics is completely described by HQET, which is given 
by the 1/rriQ expansion of the heavy quark QCD Lagrangian: 
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L = Q{iip-m Q )Q 

oo / y \ n 

=^+^X^) c " =Ca+c - (15) 

with 

Co = h v iv ■ Dh v , (16) 

and 

C n = h v (t P ± )(-iv ■ D) n -\i p ± )h v , (17) 

where Q(x) = e~ im ® v ' x [h v (x) + H v (x)] and tfh v (x) = h v (x), $H v (x) = —H v (x), and = 
- v^v ■ D. 

The heavy meson bound state, Eq. @, in the heavy quark limit is then determined 
by the leading Lagrangian Co in the above 1/rriQ expansion of QCD. Quantizing Cq on the 
light-front by letting h v = h v+ + h v _ (h v ± = §7°7 ± /i, ) ), we can derive the following light-front 
bound state equation for the heavy meson bound state (§): 

{2A-(X + ^±^)}^(X,. ± ) 

r rl Y' ' cP-k' 

This bound state equation was also derived within the light-front HQET |8,27], where the 
variable X is given by Ay. The effective heavy-light quark interaction Vgq can be obtained 
from the leading Lagrangian of (|1^) through the use of the light-front similarity renormal- 
ization group approach. Such an effective interaction is independent of the heavy quark's 
mass and spin. In other words, VQ q (X, X', k' ± ) respects HQS. For a detailed derivation, 
see §. 

It is easy to prove that the light-front wave functions discussed above possess heavy 
quark symmetry. Using the normalization property of Eq. (0), Eq. (|1|) can be reduced to 

{2A-(X+^±^)}$^(X,4) 

= / ^^^( X '^' X ^^)^ 5Sz ^( X ^")' ( 19 ) 



where 



M ss > = X - — _L=_ = Tr\^-l)v{p,X g )v(p',X q )], (20) 



yj(v-p q + m q )(v-j/ q + m q ) 



which is independent of the mass and spin of the heavy quark. Thus, the integral kernel, 
VQq(X, X', k' ± )A4 SSz which determines the wave function $ s ' s ' 2 , is also independent of 
the mass and spin of the heavy quark. Therefore, the wavefunction <& SSz (X, k 2 ± ) must be 
degenerate for S = and S = 1. As a result, we can simply write 

$ ss *(X,k 2 ± ) = $(X,k 2 ± ) (21) 

in the heavy quark limit. Eq. (§) together with Eq. (^l|) is then the heavy meson light-front 
bound states in the heavy quark limit that obey HQS. 
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D. The covariant form and diagrammatic rules 



The above light-front heavy mesonic bound states in the heavy quark limit manifest 
boost covariance and heavy quark symmetry. In fact, it can also be rewritten in a fully 
covariant form if the light-front wave function <3?(X, n\) is a function of v -p q (see the proof 
given in Sec. IV): 

—►$(*;. p ff ), (22) 
where the antiquark q in bound states is on-mass-shell, p~ = -^f(p q ± + m^). Hence, 

v-p q = ^^l + m 2 q + X 2 ). (23) 

As we will see later, the widely used Gaussian-type wave functions have such a structure in 
the heavy quark limit. 

When the wave function possesses such a form, the heavy meson bound state in the HQS 
limit, Eq. can be expressed in a covariant form: 

\H(v,S,S z )) = h v T H q v \0), (24) 

where T H = i-f 5 for S = and — ^ for S = 1, and 

K{x) = f ^(2ir)6(2v ■ k) £ {«(«, X Q )b v (k, Xq)e~ ix k + v(v, X Q )dl(k, X Q )e^ k ) , (25) 



£ {«(j9„ A g )6 9 (^, \)e— ^ + v( Pq , \)d\{p q , X q )e^} , (26) 

here is the residual momentum of the heavy quark and p q the momentum of the light 
antiquark. Note that the energies of physical particles in the above integrals are always 
positive no matter if they are expressed in terms of the light-front coordinates or the equal- 
time coordinates. This is because we have already separated explicitly the particle and 
antiparticle components in its Fourier transformation. On the light-front, this implies that 
n ■ p q = p q > [ n is a light-like unit vector defined as n M = (0, 2, OjJ] in all the subse- 
quent calculations presented in this paper. This restriction does not cause any trouble for 
maintaining the covariance. 

The Fourier transformation of h v (x)Tjjq v (x), 

(h v T H q v )(v) =vnj dx~d 2 x ± e^™ h v (x)T H q v (x) , (27) 

leads to Eq. (|8|), where the factor v ■ n in front of the Fourier transformation is due to the 
required normalization of the bound state, Eq. (|i~3P . Now, the normalization condition of 
Eq. ( |I4"D can be written in a covariant form: 
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The left hand side of the above equation can be easily obtained in a diagrammatic way as 
shown in Fig. 1: 

Fig. 1 = / ^fj^m - ml) 

= Eq. (28) . (29) 

In general, the diagrammatic rule within the covariant model is given as follows: 
(i) The heavy meson bound state in the heavy quark limit gives a vertex (wave function) 
as follows: 



y 1 -—*{v- Pq )r H , (30) 

2 V v ■ p q + m q 



2 V V ■ Pq + m q 



**(v-p q )r H , (31) 



with a momentum conservation factor 2(27r) 3 v ■ n 5 3 (v A — k —p q ). Here we do not have the 
four-momentum conservation since the meson and valence quarks are all on-mass-shell, but 
this does not affect on the covariant calculations we will perform later. 

(ii) The internal line attached to the bound state gives an on-mass-shell propagator, 



1+ ff (for heavy quarks), (32) 
i>q — m q (f° r light antiquarks), (33) 



Pq 

where v 2 = 1 and p 2 q = m 2 . 

(hi) For each internal line attached to the bound state, sum over helicity and integrate 
the internal momentum using 

J -0^(2n)5(2v ■ k) and J ^(2n)8(p 2 q - mj) (34) 

for the heavy and light quarks, respectively, where the delta functions come from the on- 
mass-shell condition. 

(iv) For all other lines and vertices that do not attach to the bound states, the diagram- 
matic rules are the same as the Feynman rules in the conventional field theory. 

These are the basic rules for the subsequent evaluations in the covariant model. Now we 
can see that the the heavy meson states in the heavy quark limit consists of a bare heavy 
quark h v surrounded by a brown muck denoted by q v , as shown in Fig. 2. The light-front 
wave function $(v-p q ) here can be regarded as the effective structure of brown muck resulted 
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from the complicated low-energy interactions of quarks and gluons. Therefore the physical 
picture of our covariant light-front model is very clear. 

Once the light-front bound states in the heavy quark limit is constructed, all the results 
obtained from the above bound states are ensured to be consistent with HQS. Furthermore, 
from such bound states, we can further develop a systematic approach within HQET to 
calculate the I/tuq corrections in the first-principles QCD. 



E. 1/rriQ corrections 



As we have discussed in Sec. II. B, the heavy meson bound states in the infinite quark 
mass limit is dynamically determined by the leading Lagrangian £ m the 1/rriQ expansion 
of QCD. Then, the 1/mg corrections to the bound states can be easily formulated in per- 
turbation theory by taking C m as a perturbation in HQET. Consider the heavy meson state 
\H(v, S, S z )) as a unperturbated state in the light-front time x + = — oo. The physical state 
\H(P, S, S z )) at x + = is then given by 

\H(P,S,S Z )) = T+ expj - t -J^dx + J dx-d 2 x ± (jC m + £ u )^\H(v,S,S z )) , (35) 

where T + is a light-front time-ordering operator, and Lu the interaction part of the QCD 
Lagrangian for light quarks. Note that the above light-front time-ordering perturbation ex- 
pansion can be equivalently expressed in terms of the usual equal time-ordering perturbation 
expansion since the unperturbated state \H(v , S, S z )} has a covariant form. 

Formally, the 1/tuq corrections to physical observables are given as follows. Let the 
physical observable be O = Oo+O m , where Oq is the leading operator in the I/tuq expansion, 
and O m = Y,n°=i (^) On is the mQ-dependent terms in terms of I/tuq expansion. Then, 



(H(P',S',S' z )\0\H(P,S,S z )) 



\2 



= {H(v\S\S' z )\T + \o + \^- J cfix^uMOo J d A x 2 C l {x 2 ) 
+ iz!l f d 4 x 1 jC lI (x 1 )0 1 + ---\\H(v,S,S z )) 

m Q J J 

= /o + ^-/i + -"- (36) 

Since the "unperturbated" bound state has a covariant form, all the subsequent 1/rriQ 
corrections can be evaluated in a diagrammatic way, which will be illustrated in Sec. V. 

As a result, our covariant model for heavy mesons contains the bound state in the heavy 
quark limit plus a systematic approach to calculate the 1/toq corrections within HQET. 
We expect that this model can serve as a quasi-first-principles QCD description of heavy 
mesons. In the subsequent sections, we shall give some basic applications of this covariant 
model, including the calculation of the Isgur-Wise function, the decay constants and the 
heavy meson masses up to the corrections. 
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III. A CO VARIANT EVALUATION OF HEAVY MESON STRUCTURE IN 

HEAVY QUARK LIMIT 



In this section, we shall present an extremely simple evaluation of the Isgur-Wise function 
and the decay constants in the HQS limit within the covariant light-front model. 



A. Isgur-Wise functions 

In the infinite quark mass limit, the transition matrix elements of B — > D, and B — > D* 
decays are given by 

{D(v')\h c v ,Th b v \B(v)) and (D*(v',e*)\h c v ,Th b v \B(v)), (37) 

where T is a Dirac 7-matrix representing for the electroweak current. Contrary to the 
previous work, here we do not need to choose a particular Lorentz frame to calculate the 
above matrix elements. Also once we have a covariant light-front model, we can directly 
evaluate the above matrix elements without restricting to the plus component of currents. 
Using the HQS-embedded covariant bound states, Eq. the hadronic matrix elements 

of B — > D and B — ► D* decays are given by (diagrammatically, see Fig. 3) 

(D(v')\h c vf Th b v \B(v)) = Tr{ 7 5(^)r(^) 75 ^} , (38) 

(D*(v',e*)\r vl Th b v \B(v)) =ro{^(!±^)r(^) 7 B^}, (39) 
where M. is the transition matrix element for the light antiquark (brown muck): 

M=f ^(2*)5(pl-ml)V(v'-p q )*(v.p q ) m \ A (40) 

J y Z7l > \/ [V ■ Pq + «\) W ' Vq + ™, q ) 

The structure of M. dictated by Lorentz invariance has the form p[ 

M = A + B^ + Cf + D^f. (41) 

This covariant decomposition allows us to easily determine the coefficients A,B,C,D in our 
covariant model with the results: 

A= (2^ {27l)6 ^ ~ < ] li = u \ \ m « ' (42) 
J y(v ■p q + m q )(v' -pq + mq) 

d% ,„ , xf 2 2^ ■Pq)®(v-p q ) 



(2ti") 4 q ^J(vp q + m q ) (v ' ■ pq + m q ) 

lf (v + V') -p q (V- V') ■ Pq \ 

2\(l+vv') + (1-W) j' {A6) 
J y Z7l > y/{v ■p q +m q )(v' -Pq + niq) 

1 f (U + V') -p q (V- V') ■ p q \ 

2\(l + vv>) {l-v-v')!' 1 1 

D = . (45) 
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Then Eqs. (]38|) and (|39| ) can be simplified as 

(DMKr^lBCT;)) = £(v ■ ^{75(^)^^)75} , (46) 
(^Ke*)|K^|S(^))=^-^Tr{^(^)r(l±^), 75 } ! (47) 

where £(v ■ v') is the so-called Isgur-Wise function given by 
£(vv') = A - B - C + D 

lq (27r)5(p 2 -m 2 q )$*(v'-p q )<S>(vp q ) q , K == ■ (48) 



(2vr) 4 '' ' ^J( v . pq j r m? ) (V . p g _|_ m ^ 

It is easy to check that at the zero recoil v ■ v ' = 1, i.e., v' = v, we have 

«1) = / ^(^(p! ft) | a = 1, (49) 

as a result of the normalized wave function, which is model independent as required by HQS. 

Since v ■ p q = ^(^± + m 2 q + X 2 ) and v ' ■ p q = ^(/s'l + m l + X' 2 ), the integration in 
Eq. ( f48[) can be explicitly expressed by 

X'4 + X/t'l + (X + X')(m 2 + XX , ) + 2XX'm (? (l + u • u') 

; , (50) 

XX'[«i + (m, + X) 2 ] [k^ + (™, + X') 2 ] 

where X' = zX, k' ± = k±_ + X(t>j_ — j_) with z = v + /v . To simplify the result, we can 
choose v± = v ' ± = without loss of generality, then z can be related to v • v' by 

1 



z — > z± = v • v' ± J (v ■ v') 2 — 1 , z + = — , (51) 



where the + (— ) sign corresponds to v 3 greater (less) than v' 3 . It turns out that 

k\ + (m q + X){m q + zX) 

x sfW± + K + x ) 2 \ [«i + K + **) 2 ] 

= c(«) = C(1A) • (52) 

Note that v 3 greater (less) than i/ 3 corresponds the daughter meson recoiling in the negative 
(positive) z direction in the rest frame of the parent meson. In other words, after setting 
v± = v' ± = 0, the daughter meson recoiling in the positive and the negative z directions are 
the only two possible choices of Lorentz frames. The last equality of the above equation 
shows that the Isgur-Wise function thus obtained is independent of the recoiling direction, 
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namely, it is truly Lorentz invariant. The above Isgur-Wise function is an improvement to 
that given in |27 where the calculation is somewhat oversimplified. 

We also note that the above result for the Isgur-Wise function is the same as that 
extracted from the plus component of P — > P transition matrix element at time-like q 2 in the 
non-covariant LFQM However, in our previous work based on the LFQM pi |, the Isgur- 
Wise function extracted from the plus component of P — > V matrix element is different in 
its expression from that for P — > P transition. This indicates that the usual LFQM utilized 
before is not fully consistent with HQS. In other words, in the usual LFQM description, 
extracting the hadronic form factors from the plus component of hadronic matrix elements 
may lead to ambiguities due to lack of Lorentz covariance. Therefore, the LFQM commonly 
used in the literature for the calculation of various form factors may not be trustworthy. Here 
we have shown that the Isgur-Wise function has exactly the same expression for P — > P and 
P — > V transitions in a simple covariant description. Our covariant model removes possible 
ambiguities in the usual LFQM calculations. 



B. Decay constants in the heavy quark limit 

To further examine the covariant model, we shall evaluate in this section the decay 
constants in the heavy quark limit. 

The decay constants of pseudoscalar and vector mesons are defined by (Ol^jP) = ifpp^ 
and (0| V^\P*) = f v M v e^, where = qj^jsQ and = q^^Q, with q and Q the light and 
heavy quark field operators, respectively, 

q{x) = I j^r {2Tr)6{p « - m " ] £ W^' XqMPq > Xq)e ~ lx ' Pq 

+ v(p q ,X q )dl(p q ,X q )e ix ^}, (53) 
Q( x ) = e - im ^ x h v + 0(l/m Q ) , (54) 

and the Fourier transformation of h v in the momentum space is given by Eq. (|25|) . Also 
note that in Eq.(|53D, n ■ p q = p+ > 0. In the infinite quark mass limit, the decay constants 
have the expressions 

(0\qr^h v \P(v)) = iF P v» , (Ol^^l^e)) = F v e» , (55) 

so that 

F P = f P jM~ P , F v = f v ^. (56) 

HQS demands that F v = F P . 

Now, using the bound states in Sec. II, it is very simply to evaluate the above matrix 
elements (diagrammatically shown in Fig. 4): 

(0\qr l5 h v \P(v)) = -^{7^75^75^1} , (57) 
(Q\7fh v \F t (v f e)) = Tt{f^ /Mi}, (58) 
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where 



and 



Mi = I j5L(ar)*G* - "S)J=^K- A) = a. + bi y . (59) 



A >=fcf^w-^^!0=>»., (60) 

* = -^/ Sy^« - m «)^Pk " ' p " (61) 



Here N c = 3 is the number of colors. Thus, it is easily found: 

F P = 2{A X - By) = 2^N C J —^(2ir)8(p 2 q _ m 2 q )<l>(v ■ Pq )y/vp q + m q 



as expected from HQS. Note that this result is different from the non-covariant LFQM 
expression given in ]rjj where again, due to the non-covariance, the decay constant is in- 
correctly extracted from the plus component alone of the corresponding current matrix 
elements. 



IV. NUMERICAL TEST OF THE CO VARIANT LIGHT-FRONT MODEL 

In the previous sections, we have constructed a covariant light-front model and evaluated 
in a very general but tremendously simple way the Isgur-Wise function £(t> • v') and the 
decay constants Fp and Fy in the heavy quark limit. As we have pointed out, covariance 
requires that light-front wave functions be a function of v -p q . Currently, there exist several 
phenomenological light-front wave functions commonly utilized in the literature. In this 
section, we shall show why a covariant light-front wave function must be a function of v -p q , 
and see if currently used wave functions obey this covariant condition. Then we shall give 
a numerical calculation of the Isgur-Wise function and decay constants with the covariant 
wave functions. 



A. Covariant light-front wave functions 

Careful readers may wonder what does covariance mean for a bound state ? Eqs. (PH), 
(pop and (J26D are only recast in a covariant form, but the internal momenta of the quarks 
carried inside mesons are all on-mass-shell. After the integration of the light-front energy is 
performed, the bound states are apparently not covariant any more. We should emphasize 
that the validity of the covariant conditions given by Eqs. (fy]) and ( |59|) depends on the spe- 
cific form of light-front wave functions. For an arbitrary light-front wave function $(X, k\), 
whether Ai and A4i have the covariant structure shown in Eqs. (]41|) and (j59|) respectively, 
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depends on the specific form of light-front wave functions. We find that to obey Eqs. (|4lD 
and (|59|) , the following identity must be satisfied: 

Eq.(|63D can be satisfied if <3?(X, k\) is a function of X + m 3+ K ± = 2v ■ p q . This leads to the 
covariant condition Eq. (p2j). 

The condition ( |63D can be easily obtained by considering the integral involved in the 
calculation of Eqs. ( fUD and fl55|): 

/ |^ (27r)( ^ ~ ^ = avAt ' (64) 

where the right hand side is the requirement of covariance, and the constant a is given by 

a = J ^(2t)*(p5 - m > ■ PMX, 4) • (65) 



Using the light-front relative momentum, the integral (|64|) for each component of p^ gives 

/ 0j; {27r)6{p « ~ m >t*( x > K ^ = V+ I ^Jc x ^ K ^ > ^ 
/ |^ (27r) ^ - m >«±*(x, 4) = «i / 4) , (67) 

f d A p a , xr/9 9x 9 , /■ dXd 2 K ± ( v 2 ! K i+ m » 1 \x/„ 9x / x 

/ (# (2 ^ (p » - m » )p « * ( - Y - = / m4 ( tj+ + ) •<*• • (68) 



It is easy to check that the above integrals can have the covariant form of (p^ ) only if Eq. (|63] 
is satisfied. Then, Eq. (|64j) can be written in a nice covariant form, 



' l pq i27r)S(p 2 -m 2 )p^(v-p g )=v^ [ ^h-(2n)6(p 2 - m> • Pq $(v ■ p q ) . (69) 



(2tt) 4V ' q/rq v 7 (2tt 

There are several popular phenomenological light-front wave functions that have been 
employed to describe various hadronic structures in the literature. Two of them, the Bauer- 
Stech-Wirbel (BSW) wave function $ B sw( x , 4) |23] and the Gaussian-type wave function 
^g{x,k 2 ± ) ||28fl , have been widely used in the study of heavy mesons. In the heavy quark 
limit, we denote these wave functions as follows ||27|J2l||: 



<W(X,4) = ^^(^)exp{-^^^}, (70) 



iG( x,4) = 4(^) s/4, - Y + m ' + K; 



V 2X 2 
XeXP |-2^ 



4 + (X/2 - K + 4)/2X) 2 . (71) 
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Since we have used the fully boost covariant form for the light-front bound states which 
are different from the previous work, the wave functions <&(X, 4) in the present paper are 
related to &(X, 4) employed in the previous works |27]J21| as follows: 



(72) 



<$>BSW(X } = VX $BSW(X, K 

$ G (x, k 2 ± ) = Vx$ g (x, 4; 




2X 



x exp 



2u 2 



4 + (X/2 - {m 2 q + K 2 ± )/2Xf 



(73) 



There is a variant of the Gaussian-type wave function which is simply a Gaussian function 
of the light-front (boost) invariant mass Mq. Unfortunately, this wave function does not 



have an appropriate heavy quark limit ||21|| . However, by comparison with the bound state 
equations (0) and fll"8|) in Sec. II, we see that the residual light-front (boost) invariant mass 

in HQS limit is given by the form of (X + K± „ mg ) = 2v -p q . Thus, we propose here a similar 
Gaussian-type wave function in terms of the residual light-front invariant mass in the HQS 
limit as follows: 



$ M pT, 4) = AT exp I - — \X + 



2u 



4 + m q 



X 



M exp 



V-Pg 



LO 



(74) 



where M is a normalization constant. 

One can see that not all the phenomenological light-front wave functions have the co- 
variant property. We found that the Gaussian wave function and the invariant-mass wave 
function can be reexpressed as a pure function of v ■ p q . The wave function $m has already 
been in this form, and <3?g can be rewritten in terms of v ■ p q : 



$ G (X, k] 



IT ^3/4 JX 2 +m 2 q + 



4 



LO ~ 



IX \3/4 



2X 



x exp 



y/v ■ p q exp 



2lu 2 



{X/2 + {m 2 q + Kl)/2X) 2 -m\ 



2^ 



{v-Pqj 2 ~m 



(75) 



Therefore these two wave functions preserve the Lorentz covariance of Eqs. (41) and 
as we have numerically examined from Eq. (|63|). However, very surprisingly, the commonly 
used BSW wave function cannot be recast as a pure function of v-p q , and it does not satisfy 
Eqs. ( PD and flSTJp. Hence the BSW wave function breaks the Lorentz covariance. Indeed, 



we have already found in a recent work |2j] that there is some inconsistent problem by using 
the BSW wave function to calculate various transition form factors. Now we can understand 



why the BSW wave function gives such results inconsistent with HQS found in Refs. [21,22 



It is interesting to see that by demanding relativstic covariance, we can rule out certain 
types of heavy meson light-front wave functions. 
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B. Numerical calculations and discussions 

Now we present some numerical results. As we have seen the phenomenological wave- 
functions have an unknown parameter uj. This parameter is proportional to the averaged 
relative momentum between the heavy and light quarks, J (/«j_)> which corresponds to the 
so-called "Fermi motion". It is of order A and can be fixed by the decay constant fg and 
light quark masses. For example, with the input 

m 9 = 0.25 GeV , f B = 0.180 GeV, (76) 

we obtain the value of uj for Q M and $g listed in Table I. We then use these values to 
calculate the Isgur-Wise function and its slope p 2 at the zero-recoil point (see Table I and 
Fig. 5). 

Table I. Parameter uj fitted to the decay constant fg with m q = 0.25 GeV in two different 

covariant light-front wave functions. Also shown are the results for J and the slope p 2 of the 
Isgur-Wise function at the zero-recoil point. 





uj (GeV) 


sJ{kI) (GeV) 


fs (GeV) 


P 2 




0.334 


0.419 


0.180 


1.09 




0.490 


0.490 


0.180 


1.26 



It is interesting to note that, because of the loss of covariance, the BSW wave function 
fails to give a finite result for fg. This reinforces the statement that one cannot use a 
non-covariant light-front wave function in practical applications. 

From the above numerical results, we see that both the covariant wave functions $m 
and <3>g give a reasonable behavior of the Isgur-Wise function. Quantitatively, different 
phenomenological wave functions give slightly different results. To get some insight of which 
wave function may be more reliable, we have also looked at the light quark mass m q and uj 
dependence of p 2 , where p 2 has been measured experimentally. The results show that for 
both $ M and $g wave functions, the slope p 2 of the Isgur-Wise function increases with a 
decrease of the averaged relative momentum \k±\, but the response of p 2 to \k±\ is more 
sensitive for $c than that for Furthermore, the current experimental analysis gives 

p 2 ~ 1.0 [29] which seems to favor the wave function $ M • Of course, there are many other 



possible wave functions and the real solution should be determined from Eq. ([Tq ) which 
deserves further investigation. 



V. 1/Mq corrections to heavy meson masses 

In this section, we present an example, the mass splitting of pseudoscalar and vector 
mesons, to show how to calculate 1 / mq corrections in our covariant light-front model. 
In the heavy quark limit, the heavy meson masses can be written as 

M H = m Q + A. (77) 
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Here A is a constant (independent of mq) and is of order Aq C d an d y(/Cj_). In principle, 
A can be directly determined by solving the bound state equation in the heavy quark limit 
from the leading order Lagrangian, namely Eq. flTBp . This is beyond the scope of the present 
paper. 

However, our covariant model permits a systematic evaluation of various I/tuq correc- 
tions, as we have formulated in Sec. II. In the infinite quark mass limit, the pseudoscalar 
and vector mesons are degenerate. Of course, in reality they are not. Their mass difference 
comes mainly from the HQS-breaking leading 1/mq corrections. From Sec. II, these I/tuq 
corrections can be expressed as 

M H = m Q + A + -^—{\ 1 + d H \ 2 ). (78) 

In the covariant model, Ai and dn^2 can be diagrammatically depicted as shown in Fig. 6: 
As we can see, Ai is the same for both pseudoscalar and vector mesons, whereas dn^ 2 
gives rise to the mass splitting between these two states. Hence, we shall only focus on 
the calculation of dn\ 2 . From the diagrammatic rules of the covariant model, we can 
immediately write down the expressions of dn A2 as follows: 

d H X 2 = Trjr^^a^^r^Mf ) , (79) 



and 

M 2 ^ = -9 2 C f J^J ^n)S(pj q - m;)(2*)*04 - m\ 



x 



r ${V ■ p lq )®*{v ■ P2g ) 

{v-pi q + m q ){v-p 2q + m q ) 

x q^D^'Up 1 ^ - m q ) lv ,(fo q - m q )) , (80) 



where Cf is a color factor, C/ = N £ N 1 = 4/3 for N c = 3, and D^ v is the gluon propagator, 



in Feynman gauge, D^ v = 2 !f m 2 with q = p Xq — p 2q , and m g is considered as an effective 
gluon mass in the low-energy domain of QCD. 
Using the identity 



we obtain 



where 



iii I±^, = 0, (81) 



M 2 V -> -a^X 2 , (82) 



(2tt) 4 7 (2tt) 4V ' ^ Vq q,y ' yrzq q 'q 2 -m 
V ( v ' V\q + m q )(v ■ p 2q + m q ) 



x 7^{(v ■ Pi q + m q ) [(q ■ p 2q ) -{v-q)(v- p 2q ) 

- (v ■ p 2q + m q ) [{q ■ pi q ) -{v-q){v p lq )] } . (83) 
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Since 




we then have 

, \ , f — 3 S=0 for pseudoscalar mesons, / oc \ 

d H X 2 = A 2 x a 1 f + ( 85 ) 

11 b=l tor vector mesons . 

This is the most general result valid in any Lorentz frame. To carry out the numerical 
calculation, it is convenient to take the rest frame v — 0. The numerical result for the wave 
function $m with m b = 4.8 GeV, m g = gives a s = 0.29 by fitting to the experimental 
B-B* mass splitting AM BB * = M B , -M B = 0.046 GeV |HJ, and a s = 0.30 (with m c = 1.6 



GeV) to the D-D* mass splitting AM DD , = 0.142 GeV For the wave function $ G , 
we obtain a s = 0.30 and 0.31 for B — B* and D — D* mass splittings, respectively. 



VI. CONCLUSIONS AND PERSPECTIVE 

In conclusion, we have constructed in this paper a new light-front model for heavy 
mesons. The main feature of this model is the preservation of covariance in the LFQM. 
We show that the covariant light-front wave function must be subject to the constraint of 
being a function of v ■ p q . As a result, we find that the Gaussian- type light-front wave 
function satisfies this constraint but the commonly used BS W wave function does not. This 



explains the puzzle encountered in [^TJ why we cannot get the correct normalization of the 
Isgur-Wise function at zero recoil extracted from P — > V transition using the BSW wave 
function. The covariant model removes some ambiguities often occurred in the usual LFQM 
calculation due to lack of covariance, releases the restriction of light-front calculations that 
hadronic form factors can only be extracted from the matrix elements of the plus component 
of the corresponding heavy quark currents at the zero momentum transfer. Meanwhile, as 
we have shown in this paper, the covariant model allows us to perform an extremely simple 
evaluation of various heavy meson properties. The typical examples given in this paper are 
the calculations of the Isgur-Wise function, the decay constants and the heavy-meson mass 
splitting which is a 1/rriQ correction in HQET. Furthermore, we have also proposed a new 
covariant light-front wave function, Eq. ([74]), from the residual light-front (boost) invariant 
mass, which gives good numerical results. Further applications to other properties of heavy 
mesons and physical processes will be presented in subsequent papers. 

Another significance of the covariant model is that it provides a quasi-first-principles 
description of the heavy meson dynamics. Namely, although the heavy meson bound state 
in the heavy quark limit is still phenomeno logical in nature at present, it is constrained by 
HQS and HQET. In the meantime, it offers a systematic approach to calculate the I/tuq 
corrections based on the first-principles 1 /mg expansion of QCD. This resembles very much 
the situation of the QCD analysis of deep inelastic scatterings for light quark systems, in 
which the low energy dynamics (described by parton distribution functions) is determined 
phenomenologically and perturbative corrections are given in a fully first-principles way. 
Here the situation may even be better since the phenomenological part is constrained by 
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HQS and HQET, whereas the nonperturbative QCD dynamics should be much simpler in 
heavy quark limit, as has been discussed in Ref. ||. 

Besides many other applications, there are some further development for this covariant 
light-front model. One is the extension to heavy baryons. The other is to solve the bound 
state equation, ([PSD, in the heavy quark limit in the low energy HQET. Of course, once 
this equation is solved, the covariant model constructed in this paper is no longer a model 
description, it will become a fully QCD description of heavy hadrons. We shall investigate 
these problems in the near future. 
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FIGURES 

FIG. 1. A diagrammatic form of the heavy meson state normalization. 

FIG. 2. A schematic picture of a heavy meson state. 
FIG. 3. The diagram for B — > D, D* transitions in the covariant model. 

FIG. 4. The diagram for B and D meson decays in the covariant model. 

FIG. 5. The Isgur-Wise function as a function of v ■ v' calculated in the covariant model. The 
dashed line is for the usual Gaussian-type light-front wave function Eq. (|73|), and the solid line for 
the Gaussian wave function Eq. (74) in terms of the residual light-front invariant mass. 



FIG. 6. The diagrams of 1/rriQ corrections to heavy meson masses. 
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